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Abstract 

The decay rate of isotropic and homogeneous turbulence is known to be affected by the large- 



^ ' scale spectrum of the initial perturbations, associated with at least two cannonical self-preserving 

^' 

• ' solutions of the von Karman-Howarth equation: the so-called Batchelor and Saffman spectra. The 

Q 

•'-^ . effect of long-range correlations in the decay of anisotropic flows is less clear, and recently it has 

(— I , been proposed that the decay rate of rotating turbulence may be independent of the large-scale 



spectrum of the initial perturbations. We analyze numerical simulations of freely decaying rotating 



^^ ' turbulence with initial energy spectra ~ /c^ (Batchelor turbulence) and ~ k"^ (Saffman turbulence) 

^— ^ , and show that, while a self-similar decay cannot be identified for the total energy, the decay is 

fvq \ indeed affected by long-range correlations. The decay of two-dimensional and three-dimensional 



modes follows distinct power laws in each case, which are consistent with predictions derived 



I . from the anisotropic von Karman-Howarth equation, and with conservation of anisotropic integral 



quantities by the flow evolution. 



I. INTRODUCTION 

Turbulent flows subject to background rotation are an important problem in fluid me- 
chanics as several systems are affected by rotation. Rotation affects the Earth atmosphere 
and the oceans at large scales, it is crucial in many engineering flows such as in turbomachin- 
ery, and is also important in many astrophysical systems such as stellar convective regions 
and gaseous planets atmospheres. As many applications arise in nature and technology, it 
becomes important to understand detailed properties of these flows. 

In particular, a long-lasting problem in fluid mechanics is that of the decay of turbu- 
lent fluctuations. The rate of energy decay in turbulent flows is known to be sensitive to 
initial conditions, and no single universal solution to the decay appears to exist to which 
all flows converge asymptotically for infinite Reynolds number [1|. Considering isotropic 
and homogeneous turbulence at high Reynolds number, there are two well known canonical 
cases as far as the energy spectrum at large scales is concerned. Provided the cubic velocity 
correlation tensor {ui{x.)uj{x.)uk{'K + r)) decays sufficiently fast for large r, these cases cor- 
respond to the so-called Saffman spectrum for which E{k — ?■ 0) ~ Lk"^, and to the so-called 
Batchelor spectrum where E{k — > 0) ~ Ik"^. Which of these spectra is observed depends 
on the initial conditions, and other non-canonical cases may arise if the above condition on 
the decay of the correlation tensor does not hold (see, e.g., [2|). The quantities L and / 
are approximately conserved integrals for fully developed, freely decaying turbulence in each 
case (/ is often called the Loytsiansky integral), and have been used to determine the decay 
rate of the total energy (see J3|-l6| for examples). The quasi-conservation in each case is as- 
sociated with self-preserving solutions of the von Karman-Howarth equations, and can also 
be interpreted in terms of conservation of linear and angular momentum. While constancy 
of L is a consequence of linear momentum conservation |5| , constancy of / is a consequence 
of the conservation of angular momentum [71, \^. 

The presence of background rotation breaks down isotropy as a preferred direction arises 







along the axis of rotation (for detailed studies of rotating turbulence, see (SHuJ]). Flows 
subjected to rotation develop anisotropics which have been shown to impact dramatically 
on the decay of energy, and there are different and sometimes conflicting results in the 
literature regarding the decay laws followed by the energy in rotating flows. Moreover, 
theoretical arguments based on the isotropic von Karman-Howarth equations are not valid 



anymore, and must be extended to include the effect of rotation. The integrals / and L 
become tensors, and for axisymmetric flows some of the diagonal components of these tensors 
can be expected to replace I and L as approximately conserved quantities. Extensions of 
these arguments to anisotropic cases such as conducting flows with an imposed magnetic 
field, rotating, or stratified flows, have been recently derived in Refs. [Tl. Il2l. Il3|. 

The development of anisotropies in rotating flows have been reported both in experiments 
and simulations. Experiments of freely decaying rotating turbulence show a reduction of the 
energy transfer, evidenced by a decrease in the energy decay rate [ij]. An anisotropic energy 
flux with a trend towards quasi- two-dimensionalization has been recently reported 15| , and 



IS a. 



so evidenced, e.g., by an increase in the correlation lengths along the axis of rotation 



16l . llTJ I . Regarding the energy decay, experimental results for non-rotating grid turbulence 
with integral scales smaller than the size of the vessel 18| show an energy decay E ~ t~^'^ 
compatible with the theoretical result E ~ t~^/^ expected for large scale initial conditions 
with Saffman spectrum ~ /c^. In 19|, a transition from a -E ~ t~^/^ decay to a decay law 
closer to ~ t~^/^ is reported once rotation starts affecting the flow dynamics, together with 
a cyclone-anticyclone asymmetry evidenced by positive skewness which grows with rotation 
rate. The decay for bounded flows, in which initial integral scales are close to the size of the 
vessel, seems to change from E ~ t^^ to ~ t^^ when rotation is present 181 ]. 

Numerical simulations and models also show a slow down of the energy decay for many dif- 



ferent initial conditions 



11 



22 



2CH24J. together with the trend towards two-dimensionalization 



251 ] ■ and the higher cyclonic-over-anticyclonic activity 



261, 1271. The first theoretical 
and considered the decrease 



study of the decay of rotating turbulence was reported in 
of the energy transfer in the presence of rotation to explain the observed slow down in the 
decay, and to predict decay rates for the energy. The exponents reported were consistent 
with those obtained in simulations in 20|, and in the experiments in 18|. Also, a very close 
decay law E ~ t~°-^ was found in 11[ for a wide range of Rossby numbers (from very small 
to moderate values), when solving equations for inertial wave turbulence and for an Eddy 
Damped Quasi-Normal Markovian model (EDQNM3). 



25| a distinction was introduced between the energy contained in modes with 



Later, in 

kz = ^ {E2D1 corresponding to two-dimensional, "slow" or "vortical" modes) and modes with 
/c^ 7^ {E^Oi corresponding to three-dimensional, "fast", or "wave" modes). The authors 
carried a comprehensive study of the transfer of energy from 3D to 2D modes, for a wide 



range of Rossby numbers, and found a non-monotonic behavior for the decay of energy for 
large, intermediate, and small Rossby numbers. They focused on the intermediate Rossby 
range, where the maximal energy transfer between 3D and 2D modes occurs (see also 29|). 
In this range, they reported an initial growth of E2D as soon as the decay begins, that 
eventually results in a cross-over with the monotonously decaying E^ij. Reference 30| also 
reported different behavior in the time evolution of E2D and E^d, but considered only the 
decay of the total energy E = E2D + -E'sd proposing a decay that depends on the amount 
of background rotation E ~ t^"'^^\ The dependence of 7 on the rotation rate was derived 
using phenomenological arguments resulting 7 ~ Q~^. 

This variety of solutions naturally led to parametric studies of the initial condition space 
in numerical simulations. For non-helical rotating flows with initial integral scale close to 
the size of the domain 124], a decay E ~ t~^ was found in simulations, in agreement with the 
experimental results [l8|. For flows with initial integral scale sufficiently smaller than the 
domain size, the case of initial conditions with energy spectra ~ k^ was studied in detail in 



Refs. 



26 



27j, where different decay rates were reported for the two-dimensional and three- 



dimensional modes. While E^o was found to decay as in the case of isotropic turbulence, 
a slower decay was found for E2D in the intermediate- Rossby range as defined in 25|. The 
decay of two-dimensional modes can be correctly explained if only a few components of 
the Loytsiansky tensor are considered (here referred in short as I± and /|| based on the 
symmetries of the flow), and it was shown that I± remains approximately constant during 
the decay [27|. Finally, it was also found that the decay is influenced by the amount of total 
helicity present in the initial conditions, and by the initial degree of anisotropy in the flow. 
Recently, the anisotropic von Karman-Howarth equation was used to consider the decay 
of a rotating turbulent flow with initial conditions following a Saffman spectrum [1^. The 
approximately conserved integral L in the isotropic and homogeneous case, was generalized 
in the axisymmetric case to two integrals L± oc L^x = Lyy and Ly = Lzz- Assuming that 
both L± and Ly are simultaneously conserved (or I± and I\\ for the case of Batchelor initial 
spectrum), and using the empirical result that characteristic lengthscales of the flow in the 
direction parallel to the rotation axis grow linearly with time 12|, the author concludes that 
regardless of the initial spectrum being ~ A;^ or ~ A;^, both cases should decay with the same 
E ~ t~^ law, thus implying that it is not possible to set apart ~ k"^ from ~ k^ turbulence 
by measuring the energy decay exponent alone. 



In this paper we analyze whether the decay of rotating flows is affected by long-range 
correlations, studying the energy decay rate for flows with ~ k^ and ~ k"^ initial conditions. 
We first present a brief theoretical discussion of conserved integral quantities for rotating 
turbulence using a von Karman-Howarth equation which includes a Coriolis term due to 
rotation. We extend the previously derived results for Saffman turbulence [12], conside ring 



Batchelor spectrum. The results corroborate the phenomenological arguments used in 27] 
to explain the decay laws found numerically for energy spectra /c^. In the second part of the 
paper we use numerical simulations to analyze the decay of rotating Saffman and Batchelor 
turbulence in periodic domains. We focus on the regime of intermediate Rossby numbers 
(initial Ro ^0.1 25f|), and therefore we will not attempt a study of the dependence of the 
decay with the rotation rate. To the best of our knowledge, the free decay of energy for 
~ k"^ initial conditions has not yet been numerically studied for rotating turbulence, and the 
resulting decay laws seem to be relevant for experiments. We show that the decay of energy 
in two-dimensional and three-dimensional modes follows distinct power laws for the cases 
of Batchelor and Saffman turbulence. The decay laws are consistent with the predictions 
derived from the anisotropic von Karman-Howarth equation, and with the approximate 
constancy of either Lj^ or I± as the flow decays. 

II. VON KARMAN-HOWARTH EQUATIONS AND SELF-SIMILAR DECAY 

The dynamics of an incompressible fluid subjected to background rotation is described 
by the Navier-Stokes equation with the addition of the Coriolis acceleration, 

9tu -F a; X u + 217 X u = - VP + z/V^u, (1) 

together with the incompressibility condition, 

V ■ u = 0. (2) 

Here, u is the velocity field, cj = V x u is the vorticity, the centrifugal acceleration is 
absorbed in the total pressure per unit of mass V, and u is the kinematic viscosity. We 
assume uniform density and the rotation axis in the z direction so fi = Qz, with Q the 
rotation frequency. 

We briefly present the derivation of the von Karman-Howarth equation including the 
Coriolis term in order to find invariant quantities during the self-similar energy decay in 



rotating flows. To this end, we write the Navier-Stokes equation in a rotating frame in index 
notation, and evaluated at two points x and x' = x + r. 



at oxk oxi 



and 



du', 



k J 



dt dx'' ^ "' dx 



(3) 



(4) 



where u' = u(x'), V'^ denotes Laplacian with respect to the x' coordinate, and e is the 
Levi-Civita symbol. Multiplying Eq. ([3]) by u'-, Eq. (H]) by Uj, summing and averaging, we 
get the equation for the time evolution of the two-point velocity correlation tensor. 
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dx 



dx'. 



^ i^'jp) + iy (^^p') ]+H'^" (^/"O + 



dx'j 



+ V \UiUjjj ^^'m \^imn'^n^j ~r Cjmn'^n^*/ ' 



(5) 



If fi = 0, and assuming the cubic velocity correlation tensor {uiUju'^) and the pressure- 
velocity correlation {uip') decay fast enough with r when r — )■ oo, two possible integral 
conserved quantities can be obtained from Eq. ([5]). When the terms on the r.h.s. of the 
equation go to zero as (9(r~^), then it follows from the time derivative on the l.h.s. that 

L= f {u-u')dr (6) 

is conserved. This is know as the Saffman integral |5|. Multiplying Eq. ([5]) by r^, another 
possible invariant is obtained when the terms on the l.h.s. are (9(r~^), namely 



-/.^(u.uVr, 



(7) 



which is known as the Loitsyanski integral and whether it actually remains constant or not 
in decaying isotropic and homogeneous turbulence has been a matter of debate 3l|. As the 
energy spectrum for small values of k can be expanded as 



E{k ^ 0) ~ Lk^ + Ik^ + 

the integrals L or / are then expected to be quasi-invariants 
scale energy spectra of the form ~ k"^ (Saffman spectrum 



(8) 

during the decay for initial large- 
, |6| or ~ A;^ (Batchelor spectrum) 



J] respectively. Using these quasi-conserved integrals, the energy balance equation dE/dt ~ 
—E"^!"^ jl can be closed (as either l^u^ or l^u^ are constant, with / the energy-containing scale), 
and different decay laws arise for each scaling of the initial spectra. 



E{t) 



t-io/7 ifE(A;,t = 0)~A;4 



(9) 



These decay laws has been observed in numerical as well as in experimental studies (see, 

e.g., yyyyyyM) 

These isotropic integrals have been also assumed to remain constant during the decay of 
rotating turbulence with ~ k"^ and ~ k^ large-scale energy spectra, to predict energy decay 



laws 



28 



In the presence of rotation, the energy per unit of time (flux) transferred towards 
smaller scales is reduced as Q is increased. A phenomenological expression consistent with 
the observations is to assume that the flux is reduced by t he ratio of the turnover time 
r ~ //E^/2 ^o ^j^g ^ave time (2^)"^ 
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35| . resulting in 
E^ 



QP 



(10) 



Further assuming J or L are conserved results in the following decays. 



E{t) 



(11) 



r^/s iiE{k,t = 0) r^k"^ 

t~^l^ iiE{k,t = 0) ^k^ 
Although these power laws are close to the behavior found in experiments and simulations, 
as mentioned in the introduction a better agreement with observations can be obtained when 
the arguments are extended to consider the effect of anisotropy in the flow. 

In the anisotropic case, and assuming absence of long-range correlations, only the final 
term in Eq. ([5]) can contribute to the rate of change of integrals of the type 



J-ijmn / T^mTn \^i^j / ^^^j 



for the Batchelor case, or 



J^ii 



[uiu'j') di 



(12) 



(13) 



for Saffman turbulence. It is important to point out here that it is not trivial that in this 
case long-range correlations should vanish as in the isotropic case (specially for the pressure 
velocity correlations, see [ll| where these terms are computed using the Fourier transform 



of Eq. ([5]) and the Poisson equation for the pressure fluctuations). Thus, the vahdity of the 
assumptions can only be verified a posteriori from the results obtained from the numerical 
simulations. 

Assuming long-range correlations vanish sufficiently fast at large-scales, the problem is 
that of finding what components of the tensors in Eqs. (fT2|) and (TT3|) are still conserved in the 
presence of rotation. We derive here the Batchelor case (lijmn) under the same assumptions 
for the decay of velocity correlations with r as in the isotropic case, multiplying Eq. ([5]) by 
'^m.i^n and integrating to obtain 

d f 

-KT rmrn (uiu'j) dr = 

-2Vti I rmrn [euk (uku'j) + ejik {u'f,Ui)) dr. (14) 



Let's assume Qi = ^z^iz, so 



' rr„rr, iuiu'-)d: 



m' n \"'t'- 



dt, 

-2Vtz I rmTn {eizk (uku'j) + ejzk {u'kUi)) dr. (15) 

As the fiow has axisymmetry, we are interested in the diagonal components of the tensor. 
For j = i, 

^z^izk {Uku'i) = ^z (UxU'y - UyU,^) , (16) 

^zGjzk {u'k'^j) = ^z {KUy - u'yUx) , (17) 

where all terms with i = z ot j = z vanish. For m = n, introducing r_L = [rx,ry,0) and 
u± = {ux,Uy,0) we get 

■^ rl{u±u'j^)dr = 

-2 / rl_Qz[{uxu'y - Uyu'^) - {uxu'y - Uyu'^)]dr, (18) 

so the terms on the r.h.s. of the equation cancel and 

^Jrl{u^u'^)dr = 0, (19) 

or equivalently, 

I± = I r\ (u_lU^) dr = constant. (20) 



The case of Saffman turbulence was considered in [12] . The relevant components of the 
tensor Ljj for axisymmetric flows are L^x = Lyy and L^z- Similar arguments as the ones 
described above lead, for the xx and yy-components of the von Karman-Howarth equation, 
to the conservation of the integral 

L± = / (u_lu'j_) dr = constant. (21) 

Note that both /_|_ and L± resemble integrals that arise in the context of the decay of 
2D turbulence 36|. This is to be expected as the spectrum of an axisymmetric flow can be 



expanded, for the modes with k\\ = kz = 0, as 

E{k^^O,kii = 0)^ L^k^ + I^kl + ... (22) 

which is (except for differences in the dimensions) also the expansion of a 2D energy spectra. 
The resemblance can be expected as rotating flows tend to become quasi-2D, concentrating 
most of the energy in the slow modes with k\\ = 0. 

For axisymmetric turbulence and a Batchelor initial spectrum ~ /c^ at large scales, 
Eq. ([2nD leads to 

I± ~ i±_l\\u± ~ constant, (23) 

where l± and /y are characteristic (energy-containing) scales in the perpendicular and parallel 
direction respectively. Assuming 2D and 3D modes are only weakly coupled in rotating 
turbulence, from Eq. fllll]) we can write an equation for the decay of the energy in 2D modes 



dE2D E2 



(24) 



and from Eq. ([23 



dE2D ^2D ^1 



n^ 



5/2,1/2 

(25) 



dt jy2^ • 

Assuming E2D ~ t'^ it follows that |27| 

E2D ~ t-2/3 (26) 

if /|| remains constant (which is reasonable as 2D modes have no dependence on the direction 
parallel to rotation, and therefore l\\ = Iq is the vertical size of the box). 
For a ~ /c^ initial spectra, 

L_L ~ l']_l\\u]_ ~ constant, (27) 



and replacing in Eq. f l24|) leads to 

dE2D -^2D^|| 



dt L±n 

and a decay 



(28) 



E2D ~ t-'/^. (29) 

The results for the decay of energy in 2D modes in rotating turbulence can be summarized 
as follows, 

E2D(t) ~ < V , ; 3Q 

[r2/3 ifE(A;,t = 0)~fc^ 

It is worth mentioning that the decay of rotating turbulence with ~ k"^ spectra is dis- 



cussed in detail in 12|, where the constancy of L± is obtained from the Karman-Howarth 
equation. In that case, the author also considers the von Karman-Howarth equation for the 
2;2;-component of the two-point correlation tensor, and assuming cubic correlations decay 
fast enough also for that component, obtains constancy of Ly ~ ^n^l^n as well as of L±. The 
constancy of these two quantities, together with the empirical result /y = /o(l + n^t) (where 



K is a constant of order one 16|, l20|]) leads to a decay £'20 (^) ~ t~^. A similar argument for 
Batchelor ~ k"^ turbulence leads to the same decay, which would result in an impossibility 
to distinguish between ~ k"^ and ~ k'^ turbulence by measuring the energy decay exponent 
alone. 

Finally, it is interesting to point out that constancy of I± and L± can be respectively 
associated with conservation of the ^-component of angular and linear momentum (see 
[71, [l2|). In fact, the mean squared angular momentum H^ averaged over a cylinder-shaped 
volume V is 

(i/f) = 27rr/||Jx, (31) 

while the linear momentum Pz is 

(Ll) = 2tiVI\\L^. (32) 

III. NUMERICAL SIMULATIONS 

In the following we resort to numerical simulations to verify the validity of the assumptions 
and predictions discussed in the previous section. As a well resolved initial spectrum ~ k'^ or 



~ k^ is needed to observe constancy of / or L in simulations of isotropic and homogeneous 
turbulence (see, e.g., |3|), we use large-eddy simulations (LES) in periodic boxes to be able 
to have initial energy containing wavenumber ko ^ 40 with large Reynolds numbers at a 
reasonable computing cost. 

Before proceeding, we must caution the reader on how the simulations in periodic boxes 
should be interpreted. Flows in periodic boxes are sometimes considered as describing a 
homogeneous flow in an infinitely periodic domain. In this work, we will interpret instead 



a flow in a periodic box as artificially confined. Indeed, as was shown in [27|, when the 
integral scale of the flow is close to the domain size, the eddies cannot become larger and 
the energy decay laws obtained are those of confined flows. The confinement is artificial in 
the sense that the flows have no Ekman layers, as the domain has no rigid walls. Moreover, 
the confinement also results in a discrete set of wavenumbers and selects a discrete set of 
inertial waves which are normal modes of the domain, unlike homogeneous rotating flows 



in infinite domains which have continuous wavenumbers (see, e.g., |lll. l37l|). Note that due 
to the discretization of wavenumbers in the former case, the existence and number of exact 
resonances depend on the domain size |37(] (as this is aggravated when boxes with large 
aspect ratio and lower resolution in the vertical direction are used, we restrict our study to 
boxes with aspect ratio of unity). 

While these are unavoidable and inherent properties of the simulations in periodic do- 
mains, the choice of a relatively large initial energy containing wavenumber /cq ~ 40 ensures 
that as the equations are advanced in time the energy containing scale of the flow will re- 
main smaller than the box size, and that a large fraction of wavenumbers will be available 
at large-scales for the system to evolve. 

A. Equations and model 

Equations ([1]) and ([2]) are solved numerically using a dynamical subgrid-scale spectral 
LES model of rotating turbulence in which only the large scales are explicitly resolved. 
The statistical effect of unresolved scales on the scales larger than a cut-off wave number 
{kc = N/2) are modeled with an eddy-viscosity and eddy-noise that are obtained after 
solving in each time step the eddy damped quasi-normal Markovian (EDQNM) equations 



for the spectrum of the unresolved scales. Details of the LES can be found in Ref. 38 1 



and in Refs. 39|, |40|] for its extension to the rotating case. A validation of the LES against 
direct numerical simulations (DNS) for the case of freely decaying rotating turbulence can 
be found in [27|. 

The simulation domain is a 3D periodic box of length 27r with spatial resolution of 
256'^ grid points (resulting in k^ = 128). A parallelized pseudo-spectral method without 
de-aliasing is used to solve spatial derivatives, and an explicit second-order Runge-Kutta 



method is used to evolve in time 41 1. 



B. Initial conditions 



In [27[ it was observed that the energy decay is sensitive to the initial anisotropy of the 
flow (in particular, to the initial amount of energy in 2D modes). To enforce spectral and 
variance isotropy of the initial conditions, we generate at t = a random velocity field 



using the Craya- Herring decomposition 4^, |43| which associates to each wave vector k an 



orthonormal frame with axes dependant on k as 

. ,, , k X a . ,, , k X i k . „. 

a is an arbitrary but fixed unit vector which we choose to be parallel to the rotation axis, 
a = (0,0,1). This representation has been widely used to study turbulence (see, e.g.. 



44 



451). Random initial conditions are created through a superposition of harmonic modes 



with random phases 0(k), (/''(k) projected over the Craya reference system for each k, 

u(k) = e^'^i + e^-^'j. (34) 

Note that due to the geometry of the Craya-Herring decomposition there is no projection of 
u into k as the incompressibility condition stands for u(k) ± k. 

As we are interested in the decay laws followed by flows with a particular initial spectrum 
(either Batchelor's or Saffman's), we can further control the shape of the initial isotropic 
spectrum by multiplying each mode u(k) in Fourier space with < A; < /cq by an amplitude 
Ck^, where C is a constant to get initial r.m.s. velocity f/ ~ 1, and /3 is a parameter used 
to adjust the power law in E{k). For k > k^, the spectrum is continued by an exponential 
decay up to kc (the maximum wavenumber explicitly resolved in the simulations). 



C. Spectra 

In order to characterize the simulations we consider isotropic as well as anisotropic spec- 
tra. The latter is useful to analyze rotating flows in which anisotropies grow with time. 

In the simulations, the isotropic spectrum is defined by averaging in Fourier space over 
spherical shells, 

E{k,t) = ^ Y, u*(k,t)-u(k,t), (35) 

fc<|k|<fc+l 

where the star denotes complex conjugate. For the anisotropic spectrum the components of 
u(k) are integrated around the axis of rotation to obtain a spectrum that depends only on 

k±, 

E{k^) = l J2 u*(k,t)-u(k,t). (36) 

fc±<|k±l<fc±+l 

From now on we will reference these two spectra simply as E{k) and E{k^) respectively. 
More detailed information on the energy spectral distribution of anisotropic fiows can 



be obtained from the axisymmetric energy spectrum e{k\\, k±) (see, e.g., [10|, |45|, |46|). This 
spectrum is obtained after integrating the three-dimensional energy spectrum around the 
axis of rotation to obtain a spectrum that depends only on k\\ and k±, and that relates to 
E{k±) as follows: 

fc|| 
Finally, the spectrum of the initial conditions is such that the Saffman or Batchelor power 
laws ~ fc^ or ~ k'^ are imposed in E{k). The relationship between the power law in E{k) 
and of that of E{k±) can be derived as follows. Equation ( 135|) in the continuous limit can 
be expressed as 

E{k,t) = ^ f^u{Kt)edf2k, (38) 

where ^ij is the spectral tensor, i.e., the Fourier transform of the second order velocity 
correlation function $jj (mj(x, t)Mj(x + r, t)), and where f2k is the solid angle in Fourier 
space. Per virtue of Eq. (!38|) . if E{k,0) ~ k'^ then $„ ~ k'^''^. For axisymmetric flows we 
can write the energy spectrum E[k±) integrating in complex space over cylindrical shells, 

E{k±, t) = \ I *n(^±, t)k±d(Pkdk,. (39) 



TABLE I: Parameters used in the simulations: number of linear grid points A^, kinematic viscosity 
V, rotation frequency fi, Reynolds number Re, Taylor Reynolds number Re\, Rossby number Ro, 
micro-Rossby number Ro^ , and a brief description of the initial energy spectrum E{k): the power 
law followed at large-scales and the range of scales where this power law is satisfied. Values of Re, 
Rex,Ro, and Rox are given at t = 0, and at t = 20 between parentheses. 

Run N u n Re Rex Ro Rd^ Initial E{k) 

A 256 1.3 X 10-^ 1580(330) 780(100) oo(oo) 00(00) k^ (1 < k < 40) 

B 256 1.3 X 10"^ 33 1580(1090) 780(520) 0.15(0.004) 0.9(0.024) k"^ (1 < k < 40) 

C 256 1.3 X 10^"^ 1280(150) 730(60) 00(00) 00(00) k^ {1 < k < 40) 

D 256 1.3 X 10-"^ 33 1280(360) 730(180) 0.18(0.002) 1.0(0.015) k^ {1 < k < 40) 



If E{k, 0) ~ k'' then, 



E{k^, Q) = l I k'^-^k^dcPkdk, 



\f{kl + kl)^k±d(j)kdk,, (40) 

When most of the energy is in the kz = plane then, from Eq. ( l40l) 

E{k^,0) ^ ^ j kl-'d<Pk. (41) 

Even in the case in which the energy is not concentrated in the kz = plane (as in the 
isotropic initial conditions), the dependence of E{k±) will be ~ k'^'^ as it follows directly 
from Eq. (HOj) . To summarize, 

E{k) ~ r ^ E{ki_) ~ kl~\ (42) 

It is important to note that this simple relation between the isotropic E{k) and the 
anisotropic E{k±) spectra only holds for the isotropic initial conditions. However, the re- 
lation will be useful to understand the power laws observed at large scales in the following 
section. As an example, when analyzing rotating turbulence with an initial isotropic spec- 
trum E{k) ~ k'^, we will be facing the case with initial E{k±) ~ k±. Note the power laws 
E{k±) ~ k± and E{k±) ~ k^, resulting respectively for Saffman and Batchelor isotropic ini- 
tial conditions, are the ones that arise in the expansion of the axisymmetric energy spectrum 
for k\\ = and A;^ ^ in Eq. 1^. 




FIG. 1: Evolution of the isotropic energy spectrum E{k) for run A from t = to t = 20 with time 
increments At = 0.5. Note how the initial E ^ k"^ spectrum is preserved over time at the larges 
scales (smaller wave numbers). 

IV. NUMERICAL RESULTS 



We now present results stemming from four numerical simulations (see table [11). All runs 
were initialized using the random (Craya- Herring) initial conditions described in the previous 
section, with an initial isotropic spectrum E{k) ~ k'^ between < A; < 40 for runs A and B, 
and E{k) ~ fc^ between < /c < 40 for runs C and D. Runs A and C correspond to isotropic 
and homogeneous turbulence (no rotation), while runs B and D are rotating, resulting in 
homogeneous but anisotropic flows as time evolves. The purpose of simulations A and C 
is to recover well known results for isotropic ~ k"^ and ~ k'^ turbulence, and to use them 
as a starting point to analyze the more complex rotating cases. Note the decay of ~ k^ 
turbulence (corresponding to runs C and D) was studied in detail in 
here to compare with the ~ k"^ case. 



and is considered 



Parameters of the runs listed in table [Dare defined as follows: The Reynolds number is 

lU 



Re 



(43) 



where U is the r.m.s velocity and the energy-containing scale / is defined from the isotropic 
energy spectrum E{k) as 



/ = 27r 



fE{k)k"^dk 
jE{k)dk ■ 
The Reynolds number based on the Taylor scale A is 



(44) 



-Rex = — 



(45) 



10"f 
0,01 



^^^^^^^^^^^.^^^ 


: 


: 


x\ 




x\ 


r 100,00 
















: 10,00 


- X 




x\ -6/5 ; 




Vx -^ 




/\^. 


E 1,00 


$^ 






;z(t) 


x^ 




"~X-^^^^*^- 


0,10 


- -17/7 Xx 


'^ 


-10/7 X!^ 


0,01 




X-N 


x; 




10 


100 : 



0,10 



1.00 
time 



10,00 



100,00 



FIG. 2: Energy evolution for run A (solid) and C (dashed). Note the decay laws with exponents 
~ —6/5 and ~ —10/7 respectively, after an initial transient of ~ 10 turn over times. Slopes are 
shown as a reference. Inset: enstrophy decay for the same runs. The decays approach ~ —11/5 
and ~ —17/7 respectively after t ~ 10. 



where 



The Rossby number is 



A = 27r 



Ro 



jE{k)dk 
jE{k)k^dk 

U 



1/2 



2QV 



while the inicro-Rossby number is defined as 






(46) 



(47) 



(48) 



This number should be initially of order one for the scrambling effect of waves not to com- 
pletely damp the non-linear term in the Navier-Stokes equation, which would result in just 
an exponential decay of the energy 10|. 

Time-dependent parameters in table [T] are given for all runs at t = 0, and also at t = 20 
(roughly ten turnover times after the self-similar decay of energy starts). 



A. Decay of non-rotating flows 

Figure [H shows the evolution of the isotropic energy spectrum for run A, at times ranging 
from t = to t = 20. The initial ~ k^ law is approximately preserved over time at the 
smallest wavenumbers, while higher wavenumbers decay due to the effect of dissipation. This 
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FIG. 3: (a) Evolution of the integral L normalized by its initial value -^(0) for the isotropic run A 
(solid) and rotating run B (dashed), (b) Evolution of the integral / normalized by its initial value 
/(O) for the isotropic run B (solid) and rotating run D (dashed). 

behavior is consistent with the quasi-conservation of L in Eq. ([S]), for which a subsequent 
decay law E ~ t~^/^ is expected for the energy. 

Figure [2] shows the energy and enstrophy decay for runs A and C After a transient lasting 
for approximately 10 turn-over times, both runs decay with different power laws. In run A, 
the energy decays close to -E ~ t~^/^ as predicted by Eq. (jH]), and as expected from the 
shape of the large-scale spectrum in Fig. [H The enstrophy shows a faster decay close to 
Z ~ t~^^'^, consistently with phenomenological arguments which indicate that E r^ t"' leads 
to Z ~ t°~^ under the assumption of isotropy, as 



dE 

-— = -2vZit). 
dt ^ ' 



(49) 



The energy decay for run C is consistent with E ~ t"^^'"^ as expected by Eq. ([9]), while 
the enstrophy decay is close to Z ~ t^^^/^, also consistent with Eq. (H9l) . The spectrum of 
run C maintains a ~ A;^ shape for large scales (not shown). All these results are consistent 



with previous simulations and theoretical results 



m 



m. 




FIG. 4: Evolution of the isotropic energy spectrum E{k) for run B from t = to t = 20 with time 
increments At = 0.5. Note how the initial E ^ k"^ spectrum for large scales (low wave numbers) is 
preserved over time, but its overall amplitude increases. 

The prediction E ~ t~^/^ for run A assumes L remains constant during the decay. To 
verify this, and to further study whether L also remains constant in the decay of rotating 
turbulence with ~ k"^ spectrum, we show the time evolution of L for runs A and B in 
Figure [3](a). We estimated L in two different ways: A fit to the spectrum with a power law 
~ k"^ was computed for the smallest wavenumbers to obtain the multiplicative prefactor in 
Eq. (IHl) (proportional to L if the spectrum is ~ A;^), as was done in [3 |. We also computed 
the two-point correlation function from the energy spectrum using 31 1 



(u ■ u') (r) = 2 



sin(/cr) — krcos{kr) 

H/ [rC ] 7^ — r;; Gt/b, 



(50) 



'0 (kry 

and then used Eq. ([6]) to compute L. Both methods give consistent results and in the 
following we show results obtained with the last method. 

In run A, L decays to half its initial value between t = and t ~ 10, but after this time 
it only changes slowly with time and can be considered almost constant (note the energy 
between t = 10 and t = 100 changes by more than an order of magnitude). However, L 
is far from being constant in run B; the behavior of the rotating runs is discussed in more 
detail in the next subsection. 

Similar results were found for the evolution of / in simulations C and D, as shown in 
Fig. [3](b). In run C, / grows during a short transient, but stays constant during the self- 
similar decay of the energy. 

While the behavior of L and / in the simulations is compatible with the theoretical and 
phenomenological arguments discussed in Sec. [TTl the rapid growth of these quantities in the 
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FIG. 5: Energy decay for runs B (solid) and D (dashed); t~^'^ and i~^" power laws are shown 
as a reference, corresponding to phenomenological predictions for Saffman and Batchelor rotating 
turbulence based on isotropic conserved quantities. 

runs with rotation seems to question the use of isotropic integrals to derive decay laws for 
anisotropic flows. 

B. Decay of rotating flows 

The analysis for runs A and C was based on the isotropic energy spectrum and on the 
quasi-invariance of the isotropic integrals L and /. However, the growth of L (/) in run 
B (D) suggests that different arguments should be used to predict their energy decay. In 
this light, we now consider the anisotropic quantities introduced in Sec. [TTl Lj_ and I±, and 
show that further analysis can be done by studying the evolution of anisotropic spectra and 
splitting the energy decay in 2D and 3D modes. A similar method was used in 27|] to study 
~ /c^ turbulence. 

The evolution of the isotropic spectrum for run B is shown in Fig. HI for different times 
from t = to t = 20. As for run A, its initial ~ /c^ behavior for large scales is preserved 
over time. However, its amplitude at large scales increases with time, in agreement with the 
evolution of L. 

The time evolution of the total energy in run B shows that the decay is close to, but 
steeper than, E ~ t~^/^. In Figure [5] we show the energy decay for this run (we also show 
the ~ k^ case for comparison, i.e., run D). The decay deviates from the expected ~ t~^/^ 
law derived via isotropic arguments (in the same fashion, run D deviates from the expected 




FIG. 6: Evolution of the axisymmetric energy spectrum E{k±) for run B from t = to i = 20 with 
time increments At = 0.5. The initial shape £^ ~ A;_l of spectrum at large scales (low wavenumbers) 
is approximately preserved over time. 

~ t~^/^ law) and strictly speaking, a power law is hard to identify. Indeed, as was shown in 
Fig. |3l both L and I grow fast in runs B and D. As a result, arguments based on the quasi- 
conservation of L or / should be expected to fail to predict the correct energy decay. The 
reason for this failure can be associated with the fact that rotating flows are not isotropic, 
and instead become axisymmetric with a tendency towards two-dimensionalization. 

From the arguments in Sec. Ull we expect only large-scale correlations in the direction 
perpendicular to the axis of rotation to be preserved, and either L±_ or I± to remain constant 
depending on the initial conditions. Therefore, we analyze runs B and D with the aid of the 
expansion in Eq. ( l22l) and of the anisotropic energy spectrum E{k±). Following Eqs. ( HM 
and ( 130|) . we will also separate the energy into the energy of slow modes E2D, and the energy 
of fast modes E-sd- 

In Fig. |6] we show the time evolution of E{k±) in run B. Per virtue of Eq. (142|) . as 
the isotropic spectrum is E{k) ~ k"^ for small wavenumbers, the anisotropic spectrum is 
E{k±) ~ k±. As for E{k), the shape of the spectrum is preserved during the decay, and 
E{k±) ~ k±_ for small k± at all times. However, unlike the isotropic spectrum, E{k±) rapidly 
decreases its amplitude during a short transient, and then the amplitude of the spectrum at 
large scales stabilizes and varies only slowly with time. Run D shows a similar behavior in 
E{k) and E{k±), but following E{k±) ~ k^ instead. 

We saw in Sec. |TT]that constancy of L± (/_|_) may be expected for E{k) ~ k"^ (~ A;^) initial 
spectra if large-scale correlations decay fast enough, and that phenomenological analysis 




FIG. 7: (a) Ratio Lj_{t)/L±{0) as a function of time for run B. (b) Evolution of I±{t)/I±{0) in run 
D. Note that in both cases, these magnitudes remain approximately constant after a short initial 
transient. 

leads then to £'20 (^) ~ t"^/^ (£"20 ~ t~'^^^). As for L, we estimate L± using two different 
methods: By fitting the spectrum for small k±, and by using Eq. (!2T]l . The two-point lon- 
gitudinal correlation function for the axisymmetric case is estimated using Bessel functions 



and the anisotropic perpendicular spectrum (see, e.g.. 



26 



m 



K ■ O K) = 2 j E{k^)Uk^r^)dk^. 



(51) 



Both estimations give similar results and the curves discussed below are obtained from 
Eqs. f l2T]) and f lSTj) . The same procedure was used to estimate I±_ from Eq. fl20|) . 

In Fig. [7] we show the evolution of L^ and /^ normalized by their values at t = 0, for 
runs B and D respectively. Although at early times their values decrease rapidly, they 
remain afterwards approximately constant during the entire simulation (note that in the 
same simulations and at the same times, L and / increase by at least one order of magnitude). 

The evolution of the energy for run B shows interesting properties (see Fig. |8]); E2D 
initially grows until t ~ 10, when it reaches its maximum value and begins to decrease. E^d 




FIG. 8: (a) Time evolution of the energy E^d in modes with k\\ ^ (dashed), and of the energy 
E2D ill modes with ku = 0, for run B. Power laws t~^'^ and t~^'^ are indicated as a reference, 
following phenomenological arguments respectively for the decay of the 3D energy, and for the 
decay of the 2D energy based on approximate constancy of L± in Saffman turbulence, (b) Same 
for run D. The power laws t~^^'''' and t~'^'^ correspond to the phenomenological predictions for 
Batchelor turbulence. 



decreases faster than E2D and at t ^ 20 both energies are comparable. After that time, the 
system is dominated by the energy in the slow modes. Starting at t ~ 10, both energies 
show a decay compatible with power laws with different exponents. E2D{f) decays close to 
~ t~^/^, in agreement with Eq. (I5U]) and with the approximate constancy of L^, while E'^d 
decays close to ~ t~^/^, the value expected from Eq. ^ for the decay of three-dimensional 
Saffman turbulence. 

A distinct decay of E2D and E^a is also observed in run D (corresponding to ~ fc^ 
turbulence, also shown in Fig. |8]). The decay of the energy in slow modes is compatible with 
E2D{t) ~ t~^/^, as expected for Batchelor turbulence from Eq. (130|) and with the approximate 
constancy of /^ shown in Fig. [71 The 3D energy shows a decay that is close to -E3/) ~ t~^°/^. 
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FIG. 9: E2D decay as a function of time for runs B (thin) and D (thick). The expected decays 
based on constancy of L± and I± are indicated as a reference. The agreement in this case is much 
better than when assuming constancy of the isotropic integrals L and / (see Fig. [5|) . 

The decay of E2D in runs B and D is compared in Fig. [91 The two different decay laws 
followed by the energy in the two runs can be clearly identified. The result indicates that 
the decay of rotating turbulence is affected by large scale correlations, as different power 
laws can be observed for E{k) ~ k^ and ~ k^ initial spectra. 

The different power laws followed by E2D and E-^u after t ^ 10 in both runs point to a 
negligible interchange of energy between slow and fast modes at late times, as required for an 
equation like Eq. fl2lj) to hold. Runs B and D have initial values of the Rossby number which 
correspond to the intermediate Rossby range studied in [251]. In this range, maximal energy 
transfer from 3D to 2D modes takes place at early times, with E2D in some cases growing 
from energy in 3D modes as a result (see Fig. [8] (a)). After this phase, the independent 
decay of the two energies implies that these exchanges are small compared with, e.g., the 
transfer of energy from vortical motions at large scales to vortical motions at smaller scales 
in E2D, and the transfer from 3D modes to 3D modes in E^d- 



C. Energy exchange between 2D and 3D modes 



A measure of the exchange of energy between 2D and 3D modes per unit of time can be 
obtained from the flux of energy across planes in Fourier space with normal kz- The amount 
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FIG. 10: (a) Energy interchanged per unit of time S between 2D and 3D modes, (a) normalized 
by \dE2D/dt\, and (b) normalized by IdE^u/dt], for runs B (thin) and D (thick). Positive values 
of S represent flux of energy from 2D to 3D modes, while negative values indicate flux from 3D to 
2D motions. At early times S is negative, indicating energy is transferred from 3D to 2D modes. 
After t ~ 10 in run B, and f ~ 10 in run D, when the distinct power law decays are observed for 
E2D{t) and E3£){t) (Fig. [8]), S becomes negligible. 



of energy per unit of time going through any of such planes is given by 



n(fc||) 




ut- • (u • Vu). dkrdknidk 



xu,nyu,nzj 



(52) 



fc,,=0 



where the integrals in k^ and ky run over their entire range (see |47| for definitions of 
anisotropic fluxes and spectral transfer functions), and the hat denotes Fourier transform. 
The energy transferred from 2D modes to 3D modes per unit of time is then 



S = U{k\\ = 0). 



(53) 



When S is positive, it represents transfer of energy from 2D to 3D motions (i.e., energy going 
from k\\ = towards larger values of ^n, resulting in an effective source of 3D energy coming 
from 2D loses), while when negative it represents transfer from 3D to 2D motions (i.e., a 





20 40 50 80 100 120 





40 50 80 100 120 



FIG. 11: Axisymmetric energy spectrum e{k\^,k±)/ sinO at different times in run A. The circular 
contour levels indicate an isotropic energy distribution. 

source of 2D energy from 3D loses). For these sources (or loses) of energy to be negligible, 
they must be small when compared with the time derivatives dE2D/dt and dE^u/dt. When 
they are, the decay of E2D and of E^d can be considered independently. 

Figure [TO] shows the function S normalized by \dE2D/dt\ and IdEso/dtl for runs B and D. 
In both runs S is negative before t ^ 10, indicating energy goes from 3D to 2D modes, and 
explaining the growth of E2D at early times in run B. Note that in run B, \dE2D/dt\ becomes 
zero at t ^ 5, resulting in an infinite value of the ratio S/\dE2D/dt\ in Fig. [10] (a). Later, 
for t ^ 10 in run B, and for t fti 20 in run C, S/\dE2D/dt\ and S/\dE^i:,/dt\ become small, 
indicating the interchange of energy between slow and fast modes becomes negligible when 
compared with the energy decay rates, and justifying the use of separate balance equations 
to study the decay of E2D and E^r,. Indeed, to have S <^ \dE2D/dt\ and S <^ \dE^D/dt\ is 
enough to consider equations for the evolution of E2D and E^£, as was done in Sec. [Til as 
is a more precise definition of what was meant in that section by the condition of "weakly 
coupled" 2D and 3D modes. 

It is worth pointing out that the fact that S/\dE2D/dt\ and S/\dE^£,/dt\ are much smaller 
than unity does not imply that the 2D and 3D modes are completely decoupled, and only 
implies that almost no energy is interchanged between slow and fast modes during the self- 
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FIG. 12: Axisymmetric energy spectrum e(A;||, A;_L)/sin^ at different times in run B. For rotating 
flows, tlie spectral distribution of energy becomes anisotropic with more energy near the k\\ = 
plane. The cross indicates the maximum of the spectrum at the different times. 

similar decay. Triadic interactions, e.g., between two fast modes and one slow mode, can 
still occur and be relevant, as long as the slow mode is only an intermediary and does not 
receive from nor give net energy to the fast modes. Indeed, such interactions between slow 
and fast modes that preserve the 2D and 3D energy separatedly would be required to have 
dE2D/dt proportional to 1/Q in Eq. (125|) . as in their absence the decay of E2D should be 
that of 2D Navier-Stokes and independent of Q (note triadic interactions between three slow 
modes are trivially resonant and independent of Q). 

D. Energy spectral distribution 



A deeper understanding of the development of anisotropy in the flow cannot be obtained 
solely from studying the reduced spectrum E{k±) and the isotropic spectrum E{k). To 
further investigate the energy spectral distribution, we present contour plots of the axisym- 
metric energy spectrum e{k\\, k±) for runs A, B, and D in Figs. [TTl [121 and [T3] respectively. 
Note that in order to obtain circular contour levels when the spectral distribution of the 
energy is isotropic, in these figures the axisymmetric energy spectrum is divided by sin 6', 
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FIG. 13: Axisymmetric energy spectrum e{k\\,k±)/ sinO at different times in run D. As in run B, 
the spectral energy distribution is anisotropic. The cross indicates the maximum of the spectrum 
at the different times. 



where 6 = arctan(A;||/A;^). 

In the absence of rotation (Fig. [TTj) the spectrum shows an isotropic energy distribution 
evidenced by circular contour levels which maintain their shape as the flow decays. However, 
when rotation is present, the spectral distribution of the energy becomes anisotropic with 
more energy near the k\\ = axis (see Figs. [T^ and [T^ . Already at t = 10, the maximum 
of the axisymmetric energy spectrum takes place in the k\\ = axis, in agreement with the 
preferential transfer from 3D modes towards slow 2D modes observed at early times in the 
flux in Fig. [10] (as a matter of fact, the maximum of the spectrum is already located in the 
k\\ = axis for times as early as t = 5). 

After t = 10, the exchange of energy between 2D and 3D modes becomes negligible. As 
already observed in the time evolution of E2D and E^d, the energy in the 3D modes decays 
faster than the energy in 2D modes. This results in an increase of the spectral anisotropy 
as time evolves, with most of the energy near the /cy = axis at late times. At the same 
time, the peak of the spectrum slowly moves through the k\\ = axis towards smaller values 
of kj_, although even at late times the peak is sufficiently far from kj_ = 1. 

The shape of the spectrum near the k\\ = plane is of interest for many theories of 
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FIG. 14: Axysimmetric energy spectrum e{kii, k±) for run B for different values of k\\ at (a) t = 5, 
and (b) t = 10. A fcj slope is shown as a reference. 



homogeneous rotating turbulence (see, e.g., ll|, |48|]). It is important to note that here the 



domain with finite size results in a discrete set of values for k\\ so the k\\ — ?■ limit cannot be 
studied, unlike homogeneous flows in infinite domains for which k\\ is a continuum. Bearing 
in mind this limitation, we show in Figs. [H] and [13 the axisymmetric energy spectrum 
e{k\\, k±) for fixed values of k^\ (from to 5), for runs B and D at two different times (t = 5 
and t = 10). 

Figures [T4l and [T5] confirm that most of the energy is contained in modes with k\\ = (2D 
modes), and that the wavenumber k± at which the peak of the spectrum takes place slowly 
decreases with time (although the wavenumber is larger than k^ = 1). As an example, in 
run B the peak of the spectrum for k\\ = is at /c^ ~ 20 at t = 5, at k± ^ 15 aX t = 10, 
and at /c_L ~ 8 at t = 35. The spectrum for modes with k\\ = also shows a positive slope 
for wavenumbers smaller than the energy containing wavenumber, and the slope seems to 
be approximately preserved as the system evolves in time. 

The spectra e{k\\, k±) for k\\ > 1 are approximately fiat for small wavenumbers (specially 





FIG. 15: Axysimmetric energy spectra e{kH,k±) for run D for different values of kii, at (a) t = 5, 
and (b) t = 10. 

for run B), and show a slope compatible with a fcj^ power law (specially for early times, at 
t = 5, and for the spectrum corresponding to k\\ = 1). This is consistent with the scaling 
effcii, /c I ) ~ /cy /cj^ valid for k\\/k±_ <^ 1 derived from inertial wave turbulence equations in 



Ui^^) 



48|, and also observed before in numerical simulations in 



m. 



V. CONCLUSIONS 



We have analyzed numerical simulations of freely decaying rotating turbulence in periodic 
domains in the regime of intermediate Rossby number {Ro ~ 0.1), with initial energy spectra 
~ k'^ and ~ k'^. We first presented a brief theoretical discussion of conserved integral 
quantities for rotating turbulence using a von Karman-Howarth equation which includes the 
Coriolis term due to rotation, extending previously derived results for Saffman turbulence 



12| . considering also the Batchelor spectrum. Assuming 2D and 3D modes are only weakly 



coupled in rotating turbulence, constancy of the anisotropic integral quantities L^ or J^ was 



used to derive phenomenologically the energy decay rate expected for energy in 2D modes. 

In simulations without rotation, we recovered well known results for isotropic ~ /c^ and 
~ k^ turbulence, and also verified the approximate constancy of the isotropic integrals 
L and /. On the other hand, in the simulations with rotation a self-similar decay was 
difficult to identify for the total energy, while L and I were found to grow rapidly during the 
entire decay. However, the decay of energy in two-dimensional and three-dimensional modes 
followed distinct power laws, in agreement with the phenomenological predictions and with 
constancy of L± and Ij_, which were also observed to remain approximately constant in the 
simulations. 

The separate power laws followed by E2D{t) and E^oit) in decaying rotating turbulence 
require the interchange of energy between slow and fast modes to be negligible. We verified 
that this was the case by studying the fiux of energy across planes with normal kz in Fourier 



rom 3D to 2D modes, as 
but that later the fiux 



space. The results showed that energy is initially transferred 

expected for rotating turbulence at intermediate Rossby number 

of energy between 3D and 2D modes becomes negligible. Contour levels of the axisymmetric 

energy spectral distribution show that in presence of rotation the energy accumulates rapidly 

near the k\\ = plane, and that at later times the anisotropy increases as a result of the 

different decay rates of the E2D and E-^£, components of the total energy. 

Altogether, the results show that rotating turbulence is indeed affected by initial large 
scale correlations in the flow, as different decay laws arise for ~ A;^ and ~ k^ initial spectra. 
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